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Abstract 

BHAGEN-IPH is a FORTRAN program providing fast Monte Carlo event generation 
of the process e^e~ ^ e+e~7, within electroweak theory, for both unpolarized beams 
and also for the longitudinally polarized electron beam. The program is designed for final 
leptons outside a small cone around the initial leptons direction and has a new algorithm 
allowing also for a fast generation of non coUinear initial and final emission, as well as for 
asymmetric and different angular cuts for final leptons. 
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PROGRAM SUMMARY 



Title of program: BHAGEN-IPH 

Program obtainable from: the authors on request 

(E-mail: caffo@bologna.infn.it, czyz@usctouxl.cto.us.edu.pl ). 

Computer/ Operating system: any supporting FORTRAN 77. 
Programming language used: FORTRAN 77. 
Memory required to execute with typical data: about 1.6 MB 
No. of bits in a word: 32 

No. of lines in distributed program: 3607 hues. 

Subprograms used: RANLUX [9] (included in the source code). 

Keywords: radiative Bhabha scattering, polarized beam, longitudinal polarization, photon 
emission, event generator. 

Nature of physical problem: The process is measured in e+e" experiments with e~ beam 
longitudinally polarized or not, and is of interest for tests of the Standard Model (QED at 
low energies), as a background estimation for other processes, and as a radiative correction 
contribution to small angle Bhabha scattering for luminosity monitoring. 
Method of solution: a Monte Carlo event generation with importance sampling method. 
Typical running time: to generate 10^ unweighted events are requested 10-20 sec for typical 
event selection on an AlphaVAX DEC 3000 M700 with Open VMS 6.2. 
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LONG WRITE-UP 



1. Introduction. 

The cross section for e+e" i— e"'"e~7 in QED has been calculated with different 
approximations for various purposes, producing distributions to which we refer in [1]. At 
high energies (TRISTAN, LEP, SLC) it is necessary to include the weak-boson Z° exchange 
and with a longitudinally polarized electron beam (SLC) the polarization effects. For the 
unpolarized case a relatively simple expression for the square of the matrix element is 
obtained in [2,3], for the longitudinally polarized electron beam in [4] is given a reasonably 
compact expression for the square of the matrix element, and in [1] that expression is 
improved with some other terms relevant on the Z boson peak and with all the relevant 
mass-corrections for the configurations in which the final fermions have angles with the 
initial direction larger than, say, 1 mrad. For extremely forward final fermions the mass- 
corrections reported there are not sufficient. 

Because of the experimental interest in this process we have implemented it within a 
fast Monte Carlo event generator. The main target is to obtain a fast generation proce- 
dure for the mentioned radiative Bhabha scattering configurations, where final leptons are 
outside a small cone around the initial direction (of a semiopening angle of 1 mrad) and 
allowing also for the generation of a non coUinear initial or final photon without loosing the 
efficiency. At variance from the existing Monte Carlo event generators for Bhabha scat- 
tering [5,6,7], which include the radiative Bhabha scattering as a correction, we keep it as 
a distinct process and the developed method of generation allows to choose asymmetrical 
and different cuts for electron and positron angular variables, notably in t-channel. 

We do not report here the squared matrix element for the cross-section relevant in the 
region of applicability of the program, which can be found in literature for the unpolarized 
case in [1,3,4] and for the longitudinally polarized case in [1]. Here we present the details 
of the structure of the program BHAGEN-IPH, while some results and the numerical 
comparisons of BHAGEN-IPH with other existing programs (for the unpolarized case and 
special configurations) are presented in [1] . 

2. Monte Carlo Algorithm. 

The importance sampling method [8] is used as basic Monte Carlo algorithm. 
We denote the differential cross-section as 

da = 11 dQidn^dw . (2.1) 

where dQi = d(j)idcos9i {dQ^ = d(j)jd cos 9^) is the final electron (photon) solid angle, u> 
is the photon energy and S is the distribution. 

Moreover we use the following notation: dQ2 = d(f)2dcos92 is the final positron solid 
angle; Ei{E2) is final electron (positron) energy; is the initial electron and positron 
energy; rrie is the electron mass; s = 4i?^; a is fine-structure constant. We work in the CM 
frame of initial electron and positron with the z-axis chosen along initial electron direction 
and x-2;-plane given by momenta of initial and final electron. The last choice is convenient 
due to the fact that the matrix element for the unpolarized cross section, as well as the 
matrix element for the cross section with a longitudinally polarized beam, depends only 
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on one azimuthal angle, while the second one is evenly distributed. As a consequence one 
can actually generate in the chosen frame only four variables and the evenly distributed 
(in the laboratory frame) azimuthal angle of the final electron can be generated at the end 
with the subsequent rotation of all momenta. 

The distribution S is approximated by T,a-i which consists of 10 parts, in the following 
called channels, related to all possible combinations of enhancements due to Z° resonance, 
coUinear and soft photon emission and t-channel photon exchange 



S ~ = Si + E3 + Se + E7 + Eio + i^-^j (S2 + E4 + E5 + Eg + Eg) . (2.2) 

With the proper choice of the channel approximants E^, z = 1, 10, described below, 
all the relevant peaks in the distribution E can be reproduced. As each peak can be 
expressed in a more simple way using different set of angular variables, the cross-section 
(and its channel contributions exact dai and approximate daf) is rewritten as 

10 10 ^ 

da = ^dai yT'^^^ 

E (2-3) 
= —^du;\T,idQidQi^ + T,2di^2dO,2'y 

+ (E3 + Eg + E7 + Eio) dQidn^ + (E4 + E5 + Eg + Eg) dn2dn^] , 



where dQij {dQ2y) is the photon solid angle in the frame where 2;-axis is chosen along final 
electron (positron). The ratio w = (E/E^), called weight, is a relatively fiat function. 
This fact and the absorbtion of all the peaks in Ei, Eio through an appropriate change 
of variables, allow for a fast generation. Besides of the most relevant peaks, all the other 
angular and energy distributions are absorbed, whenever possible in a simple way. 

In the following are given explicit expressions for the approximants Ej, z = 1, 10, 
which define probability density in each channel. The integration for each of them is 
examined and the proper form for daf is obtained. The form integrated over one of the 
azimuthal angles is used, due to its plain distribution as previously discussed. 

We select channels according to: 

1) s-channel with final lepton emission (Channels 1-2); 

2) s-channel with initial lepton emission (Channels 3-8); 

3) t-channel with photon emission (Channels 9 and 10). 

For all the channel approximants we succeeded to obtain analytical integration (fast 
generation) and analytically invertible variables. 



1) s-channel with final lepton emission (Channels 1-2). 
- Channel 1 - final electron emission. 

To approximate well the peak in the cross-section coming from the s-channel contri- 
bution (from both photon and Z° boson exchange), when the photon is almost coUinear 
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to the final electron, we choose 

271- 



(2.4) 

Q;3 I Eh - - 

- — F{s) (l +cos^^i) —E— -j^ "^dcosOidni^dw , 

47rs ^a + 2sin2%i uj{Ei,-u) 



where 



F{s) = 2 + 44 '^'^^^'^ + 2 [(4 + 4)^ - 2P,cyCA{4 + 4)] ^ , (2.5) 

D(s) = (s-M|)2 + M|r| , (2.6) 

-1 -l + 4sin^ 9w ,^ 

4 sm t'w^ cos dw 4 sm cos dw 
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mf 



and 

o = „ 
2^6 

is Weinberg angle, Mz (F^) is Z° boson mass (width) and is the electron 
longitudinal polarization. 

Moreover the expression in (2.4) mimics the electron angular (^i) and the photon 
energetic (a;) behaviours of this part of the cross-section. The photon energy dependence 
is well approximated at both ends of its spectrum. 

By performing the change of variables 

y(^i^) = -ln(^a + 2sin2^^ , (2.10) 

^(^1) = cos^i + Jcos^^i , (2.11) 
the expression in (2.4) becomes 

do^ = — F{s)d(t)i-ydx{u;)dy{9i^)dz{9i) . (2.12) 

It is now easy to obtain analytically the integral of (2.12) over the allowed region of the 
variables, which is a key step in the importance sampling method for the choice of the 
proper channel for generation. 
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Also the relations (2.9)-(2.11) are analytically invertible: 



2 exp(x) 
exp(a;) + ■\/exp(2x) + 4 



sm 



0086*1(2;) = 



di-yiy) _ cxp(-?/) - a 



3z 



1 + 



1+1^2 



+ 



(2.13) 

(2.14) 
(2.15) 



1 + 1^2 



What exposed provides a fast generation in this channel. 
- Channel 2 - final positron emission. 

This channel approximates well the peak in the cross-section which appears in the s- 
channel contribution (from both photon and Z° boson exchange diagrams), when photon 
is almost parallel to the final positron. The situation is analogous to the channel 1, with 
the electron variables substituted by the positron ones. We use both channels to allow for 
asymmetrical and different cuts for electron and positron. 

The formulae are identical to the previous ones (with the change of subscript 1 into 
2) and the approximate expression for channel 2 is 



2w 



d(T2 = I d(l)2^2dcos92dQ,2'yduj 

(2.16) 





F{s)d4>2'ydx{(jj)dy{92'y)dz{92) ■ 



Atts 

The generation is done for convenience in the frame were x-2;-plane is given by initial 
electron and final positron momenta so after generation the momenta are rotated into our 
frame of reference (a;-2;-plane given by initial and final electron). 



2) s-channel with initial lepton emission (Channels 3-8). 

The next six channels are used to approximate the initial emission spectrum. The 
situation here is slightly complicated due to the presence of Z° boson resonance. We 
choose to use more channels and simpler formulae, instead of a numerical approximant, 
avoiding accuracy adjustments and getting faster procedure. 

Channels 5-8 are used only when the total energy is above Z°-boson mass, due to the 
fact that in the photon spectrum there is a peak, connected with the Z^-resonance, which 
is located in the hard part of the spectrum. A posteriori is found that the efficiency of 
the generator improves using these channels only when d) = (s — M|)/(4£^5) > 1 GeV. 
Covering the phase space with additional four channels allows us to get a proper shape 
of all angular distributions, avoiding the difficulty of matching more than one peak with 
the same variable. Furthermore all the variables used are analytically invertible, so the 
generation is fast. 

- Channel 3 - initial electron emission. 



6 



In this channel the influence of the Z°-boson resonance on the photon energy distri- 
bution is not taken into account. The approximant is therefore very similar to the ap- 
proximants used in channels 1 and 2 and approximates well the peak in the cross-section, 
which appears in the s-channel when the photon is collinear to the initial electron 

27r 

da^ = J dcpi'EsdcosOidn.ydLij 

(2.17) 

= — F(s) (1 + cos^ ^i) — ^ ^dcos didVt^duj , 

47rs ^a + 2sin2^ u{E^ - u) ^ 

and the peaks are absorbed in a way analogous to channel 1, with the obvious substitution 

^l/y > d-y. 

- Channel 4 - initial positron emission. 

Similarly to channel 3, the peak is approximated in the s-channel contribution, when 
photon is collinear to the initial positron, leaving the resonance as a spectator 

da^ = J d(l)2^4:d cos 92dVt^du3 

(2.18) 

3 1 JTj i£ 

= — Fis) (1 + COS^ 6*2) a- -7T=, cos e2dn^dw . 

47rs ^a + 2cos2^ u{Eb-u) 

Note however the change in the photon angular dependance, due to the fact that the 2;-axis 
is chosen along initial electron direction. Now for peak absorption and photon polar angle 
generation, instead of the variable y{9'y) in (2.10), we use the variable y{9'y) and its inverse 

2 _ 2 ^jiy) _ exp(y) - a 



y{e^) = In (^a + 2 cos^ j , cos^ = . (2.19) 

Generation is done in a frame were a;-2;-plane is given by initial electron and final positron 
momenta so again at the end we perform the rotation into the proper frame. 
- Channel 5 - initial electron emission, resonance and 62 angular distribution. 

The influence of the Z^ resonance on the photon energy is accounted for, in conjunction 
with the positron angular distribution and the photon emission collinear to the initial 
electron. As mentioned this channel is complementary to the channel 3 and of relevance 
for total energy well above the mass of the Z^. 



2tt 

da^ — I d(l)2'^5d cos 9 2dQ J du> 

o? - , . .9 1 1 16^;^ 





3 

- — Cz (1 COS 6*2)^ ^-3- - T-TT NO ^ .^2t^2 dcos92dn^dw 

47rs ^ ' a + 2sin2^ ^ 16(a; - a;)^ + M|r| 

(2.20) 
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where we introduced the notations 



± /'2, 2\2,.22 



s-Ml 



(2.21) 



With the change of variable in (2.10), with the substitution — > and the following 
ones 



5(^2) = ^(1+ cos ^2)' , 



r{u) 



■r—— arctan 



4(cD - uj)Eb 
MzVz 



MzVzCb 

the following easily integrable expression can be obtained 

>,3 



dai 



a 



c^d(f)jdy{9j)dz{d2)dr{(jj) 



These relations have the inverted expressions 

cos 02 (z) = — 1 , 



(2.22) 
(2.23) 



(2.24) 



(2.25) 



. MzTz , 
oj{r) =oj- tan 



uMzTz_ 
4EI 



(2.26) 



allowing for simple generation procedure. 

- Channel 6 - initial electron emission, resonance and 9i angular distribution. 

The influence of the resonance on the photon energy is accounted for this time in 
conjunction with the electron angular distribution and with the photon emission coUinear 
to the initial electron. Together with channels 3 and 5 gives a complete approximation of 
the s-channel cross-section with photon emission coUinear to the initial electron, when the 
total energy is above the Z^ mass. 



2n 



dai 



-I 



d(f)i'EQdcos 9idQ,^dw 



a 



c+ (1 + COS OiY 



16EI 



a + 2sm^'-^ ^ 16((I; - a;)2 + M|r| 



dcos OidVt^du) 



a 



C2;d(f)-ydy{0-y)dz{0i)dr{ijo) . 



Aits 



(2.27) 

- Channel 7 - initial positron emission, resonance and 6*1 angular distribution. 

Now the photon emission coUinear to the initial positron is considered, accounting 
for the influence of the Z^ resonance on the photon energy and in conjunction with the 
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electron angular distribution, of relevance for total energy above the Z° mass 



da^ = J d(l)iT,7dcos9idfljd(jj 





- — 07(1 — cos6'i) s- — -Tzrrz ^ ^^^^ dQ,o^d\d^-^du} 

47rs a + 2cos2% ^ 16(a; - a;)^ + M|r| ^ 

c'^d(f).ydy{9y)dz{9i)dr{u)) , 



«3 



(2.28) 



47rs 

where the variables in (2.19), (2.23) are used with the new variable 

5(^1) = -^(1- cos ^1)^ , (2.29) 

again invertible into 

cos 01 (z) = - + 1 . (2.30) 

- Channel 8 - initial positron emission, resonance and 6*2 angular distribution. 

Again the photon emission is coUinear to the initial positron, accounting for the in- 
fluence of the resonance on the photon energy and in conjunction with the positron 
angular distribution. This channel with channels 4 and 7 gives a complete approximation 
of the s-channel cross-section, when the photon is emitted coUinear to the initial positron 
and the total energy is above the Z^ mass. 

27r 

dag— d(j)2T,8d cos 92dQ^duj 







- — c^fl — cos6'2) — — ; r;^ — tt—tt d cos 9 2dO,^dL0 

Atts ^' a + 2cos2% Co 16(0; - a;)2 + M|r| ^ 
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d(i)^ dy{9^)dz{92) dr {to) 



Atts 



(2.31) 

Here the variables in (2.19), (2.23) and (2.29) are used with the substitution 9i — > 6*2. 

3) t-channel with photon emission (Channels 9 and 10). 

Even if for s-channel generation we do not follow any existing solution, the generation 
there is relatively simple. The complexity of the overlapping peaks in t-channel makes the 
problem much more complicated and the only satisfactory solution to this problem existing 
up to now [5-7] , is not efficient enough, when one is interested in generating photons outside 
small angular regions around final and initial lepton directions, which is of experimental 
relevance, when the photons are detected separately from leptons. We present the method 
we use for this case in more detail. We split the approximant of that part of the cross- 
section into two parts, which are naturally separated in the exact form of the cross-section 
and related to the emission from the electron or the positron line. 
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- Channel 9 - electron emission. 

This part approximates the peaks coming from t-channel emission from electron line, 
when photon is mostly collinear to initial or final electron. 



da 



27r 

g — J d(f)2^9d COS 62dQ^(kij 



. Ek-^ 1 



where 



—A 

TTS oj{Ei, — a;) (1 + cos ^2) a + 1 + cos ^27 a + 1 — cos 9, 



cos ^27 = cos 02 cos 9-y + sin 62 sin 9^ cos 0^ . 



-d cos 92dQydu! , 



(2.32) 
(2.33) 



The factor A will be defined later when its origin will become clear. 

With the change of variables (2.9) the approximate expression (2.32) becomes 



1 



1 



1 



^ _ ct ^ 

^ TTS 1 + COS 6*2 a + 1 + cos 6*2-7 a + 1 — cos 9^ 



dcos92dn^dx{Lo) . (2.34) 



'27 I J- ^v^oiy^ 

Splitting the angular range of into two parts [0, tt) and [tt, 27r) the variable 

2 



V{(f)j) 



arctan 



Vc2 - 62 tan ^ 
c + 6 



(2.35) 



Vc2 - 62 

can be used to absorb one of the peaks. Here the following notation is used 

c — a + 1 — cos 6*2 cos 9^ , b — — sin ^2 sin 9^ . (2.36) 
The region [0,7r) is mapped into v G [0, VjTjaa;), while [tt, 27r) into v G [ — v^naxi 

0), where 

(2.37) 



TT 



"■inu.v /— 75 =-75- • 

Vc^ — 0^ 

All that allows to generate v e [—Vmaxi'^max) with the inverse relations 



and 



(f)j = 2 arctan 



0^ = 27r + 2 arctan 



c+b 
Vc^ - b^ 



tan 



c + b 



tan 



'y/c^ — 62 



Vc2 - 62 



/or V > 



for V < . 



Vc2 - 62 

To generate in the usual interval the variable p{(f)'y) G (0, 1] is introduced 



(2.38) 



(2.39) 



P(07) 



(2.40) 
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so the approximate channel 9 contribution to the cross-section becomes 
daa = A 



^ s l + cos6'2 a + 1 — cos^-y (2 Al') 

d cos cos 9^dp{4>^)dx{Lj) 



y^a^ + 2o(l + cos 6*2 cos6'-y) + (cos 6*2 + cos6'-y)^ 

To absorb the remaining peaks and to ehminate the square root in the denominator the 
following change of variables is done 

zi (Oj) = a + l-cosd^ , (2.42) 

and 

^2(^2) = l + cos^2 , (2.43) 

so to have 



4 2a^ , 1 1 



daa = A 



dz2{02)dzi{6^)dp{(j)^)dx{uj) 



Vzf{e^) + zl{e2) - 2(a + I)z2{e2)z^{e^) - 2^2(^2) + 2(a + 1)2^2(^2) 
With the further change of variables 



D + z^{e^) + ^J {D + zi{d^)f + AP 



whose inverse is 



zi{e^)^^(^-^ -D 



The last change of variables is 



(2.44) 



t{e^, 62) = '-^^ , (2.45) 



(2.46) 



with 

L» = -(0 + 1)^2(^2) , P=^a(2 + a)^2(^2)(2-^2(^2)) , (2.47) 
the approximate expression becomes 

dai = -'^A ^^2_\t_^ dz2{e2)dt{e^,e2)dp{(t)^)dx{u) . (2.48) 



Ti{e^,e2)= y , (2.49) 



and 

^2(^2) = , (2.50) 
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where 



and 



In 



2^/Pt - D - 



2^/Pt -D + VA 



A = zi{e2)+2a{2 + a)z2{e2) , 
so that the minimum and maximum values of {m stands for min or max) are 



1 



In 



2z 



Im 



Zim - (a + 1)22(6'2) + V4m " 2-Slm(a + 1)^2 (^2) + A 



(^zim + ^/A + ^^2^ - 2zim{a + 1)^2(^2) + a)' 

and the approximate expression is now 

2a^ 



(2.51) 
(2.52) 

(2.53) 



daa 



-A {Vmax - rmin)z2dr2{02)dri{d^, 92)dp{(j)^)dx{u}) 



(2.54) 



All the variables used are analytically invertible 

1 



cos 02 = —1 + 



r-2 



COS 9^ = a + 1 — 



2Ae' 



(2.55) 



{l - e^^) VA (^e'-^ + 1) + - e*"^) 
A convenient choice for the factor A is then 

400 

— 

{rmax - rmin)z2{02) 

so to have at last the very simple approximate expression 



da^ = 



800a3 



dr2{92)dri{9j,92)dp{(j)j)dx{uj) . 



(2.56) 



(2.57) 



The choice for ^ 1 in (2.56) implies that the approximation of the peaks in t-channel 
is not perfect, but nevertheless is good enough (logarithmic in ^2(6*2)) to allow a fast 
generation. Moreover it is chosen to fix properly the relative strength respect to the other 
channels subgenerators. 
- Channel 10 - positron emission. 

This part is specular to the channel 9 and approximates the peaks coming from t- 
channel emission of a photon mostly coUinear to initial or final positron 

2n 



da 



10 



d(f)i'£iodcos 9\dVL.yduj 
1 



1 



1 



TTs oj{Ei, — u)) 1 — COS 6*1 a + 1 + cos 9i-y a + 1 + cos 9-^ 
SOOct"^ 

dr2{9i)dri{9-y,9i)dp{(f)j)dx{u!) , 



dcos 9idQ,^d(jj 



(2.58) 
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where now 

Z2{9i) = 1 - cos^i , (2.59) 

and 

zi (6*^) = a + 1 + cos 6*^ , (2.60) 

with 

c = a + 1 — cos 9i cos Oj , (2-61) 

and 

b = -sin^isin^^ , (2.62) 

while all the other expressions remain valid. A has the same functional form of A, but as 
Z2, zi, c and b are now different, it is different from A. 

Having so given the separate description of various channels, obtaining a good ap- 
proximant of the cross-section, the events can be generated in an easy way according to 
the approximated forms in each channel. 

The pseudo random number generator used to produce the equidistributed random 
numbers between [0,1) is the CERNLIB generator (V115) RANLUX [9], which source code 
is added into our program. 

The algorithm used to select the channel for generation is the following: the i-th 
channel is chosen randomly with the probability 

Pi = ^^iS^ , (2.63) 

where 

Vi= [ daf , (2.64) 

and Vi is the allowed region of the variables appearing in daf. 

A random generation of 1000 events is done at the beginning and the maximum- weight 
ws is taken to be 2.5 times the largest of the generated weights w — As the w function 
is relatively flat the method works well and in all tests an event with the weight w > ws 
was never generated. Nevertheless this eventuality is checked automatically in the program 
and in case events with w > ws appear the cross-section corresponding to them is given 
in the output. However the cross-section and its error calculated with the weighted event 
sample is not affected by the appearance of the events with weight w > ws- 

The cross-section a integrated over the allowed phase space is calculated as a sum of 
all channels contributions, according to the stratification method described in [8], 

10 10 „ 10 

^ = = E / = E / ^^^^ ' (2.65) 

i=l i=l i=l ^ 

and its statistical error Act is calculated as 



Ad 



\ 



10 



J](Aa02 , (2.66) 



i=l 
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where Acr, is the statistical error of the z-th channel contribution cr^. 

A more complicated event selection, then that allowed by choosing the range of the 
generation variables, can be implemented by simply setting the weight of the generated 
event to zero in case it is not allowed by the requested cuts (this can be easily done by the 
user in the subroutine TRIGGER). 

Finally we remark that a big effort was devoted to assure that the accuracy of the 
calculations is not lost due to possible cancellations in formulae. 

3. Routines in the Generator. 
BHAGEN-IPH 

It is the main program and calls the subroutines PARAMET, VOLUME, MAXEST, 
GENER and CROSS. It calculates and writes the total cross section for both weighted 
and unweighted events sample. 
PARAMET (il, 12, 13) 

The argument il is the logical number of the input file, from which all relevant param- 
eters are read (see input file description in the Users Guide section); 12 is the logical number 
of the first output file (BHAGEN-IPH. RES), where parameters and cross-section values 
are reported; 13 is the logical number of the second output file (BHAGEN-IPH. EVN), 
where generated events are listed (if requested). The values, which are constant during all 
the run, are calculated (such as coupling constants, boundaries of the integration regions 
etc.) and the random number generator is initialized. 
VOLUME 

The values of Vi and Pi {i = 1, 10) are calculated (see algorithm description). 
MAXEST 

It generates 1000 events in the same way as in GENER and on its basis estimates the 
maximal weight ws- 
GENER 

It chooses randomly (RANDOMNl) with the probability Pi a given channel and 
calls a proper generation procedure in each channel (GENERl, GENERIO). It writes 
the generated events (kinematics of the outgoing particles), if requested. It stops the 
generation, when the unweighted event sample is equal to the requested number. 
CROSS 

It calculates the integrated cross-sections for each channel, for both weighted and 
unweighted event sample. 
GENERl, GENERIO 

They are master programs for each separate channel. They generate (RANDOMN) 
the kinematical variables, which are chosen as independent (see algorithm description). 
The subroutines GENER3, GENERIO caU the subroutines PHASESPl or PHASESP2 
to calculate the rest of the necessary kinematics (energies and angles of the final particles) 
and reject the event if it is generated outside the allowed phase space. In GENERl 
and GENER2 the same is done inside the subroutines, as the calculations are different 
from other channels. All the subroutines call FRAME, which generate (RANDOMNl) 
the azimuthal angle of the final electron and rotates the event (up to this call the event 
is given in a frame where the a;-2;-plane is given by initial and final electron momenta). 
They also call TRIGGER, which is the event selection subroutine to be supplied by the 
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user, presently with no action. All the subroutines call FRATI03, where the weights are 
calculated and summed and where the hit or miss method is applied, selecting unweighted 
event sample. 
RANDOMNl 

It gives one random number evenly distributed on (0,1] used to choose the channel or 
to generate the azimuthal angle of the final electron. 
RANDOMN 

It gives a sequence of 5 random numbers evenly distributed on (0,1]: 4 to generate 
the event variables and 1 for the selection of the unweighted event. 
PHASESPl and PHASESP2 

They calculates energies and angles of the final particles, given 0i, ^i, 0^, 6'-y and uj 
(PHASESPl), or 02,^2,07,^7 and to (PHASESP2). 
FRAME 

The azimuthal angle of the electron is generated randomly (RANDOMNl) and the 
event (momenta of the final particles) is rotated around ^-axis by this angle. This allows 
the use of an arbitrary frame with 2;-axis along the initial electron. Before that rotation 
the a;-2;-plane is the plain given by initial and final electron momenta. 
TRIGGER 

This is an event selection subroutine, to be supplied by the user, allowing for rejection 
(k=0 has to be set inside the subroutine) or acceptance (k=l has to be set inside the 
subroutine) of the generated event. The final particles energies and angles are stored in 
the "ev" vector. They are given in CM frame of initial electron and positron, with the 
2-axis along initial electron momentum. If left as it is given in the distributed version, the 
trigger has no efi'ect on the generation. The reading of the event vector is as follow: 

ev(l) - the electron energy in GeV 

ev(2) - the electron polar angle in rad. 

ev(3) - the electron azimuthal angle in rad. 

ev(4) - the positron energy in GeV 

ev(5) - the positron polar angle in rad. 

ev(6) - the positron azimuthal angle in rad. 

ev(7) - the photon energy in GeV 

ev(8) - the photon polar angle in rad. 

ev(9) - the azimuthal angle in rad. 
FRATI03 

It calculates parameters which enter in both, the exact and the approximate cross- 
sections. It calculates the weight {w) of the event and it provides the book keeping for 
the calculation of the cross-section for the weighted sample, adding w to Wi and w"^ 
into wf for the proper channel. It checks if w is smaller of the estimated maximum- 
weight (?«5), and if not, it sums the weight w = w° into the separate sums and 
'Yliii^i)'^ allowing for the calculation of a separate cross-section of such overweighted events, 
writing also identifying parameters. The unweighted events sample can still be used, if 
the contribution to the cross-section due to the overweighted events is within the errors of 
the total cross-section. However it never happened to find an event with weight w > ws 
in the test runs, so the maximum estimation is pretty good. The hit or miss method is 
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also applied here: the randomn number r5 is generated (RANDOMN), wt = r5 ws is 

calculated and the event is accepted if w > wt- 

RFULL 

It calculates the exact differential cross-section; some factors which are the same in 
both exact and approximate form are omitted, as only the ratio of the two is used in 
calculations. 

CHANNELl, CHANNELIO 

They calculate the approximate differential cross-sections for each channel (but the 
factors which are the same in both exact and approximate form, as only the ratio of the 
two is used in calculations). 
FACT 

It is an auxihary function used by CHANNEL9 and CHANNELIO. 
4. Users Guide. 

In the version installed on VAX with VMS system the parameters are given in the 
BHAGEN-1PH.COM file (the same file where are the system commands running the pro- 
gram). In the version installed on IBM 340 (RISC/6000) with the system AIX 3.2, the 
only differences can be in the syntax of the OPEN statement (depending on the FOR- 
TRAN 77 version installed) and in the fact that the file BHAGEN-1PH.COM contains 
only parameters. The user has to supply the following data: 

- The mass of the boson and its width. 

- The mass of the W boson and the value of sin^ 6\y. Actually only one of this parameters 
is used by the program. If sin^ 0\y is set to zero the program uses the expression sin^ 9\y = 
1 — (M^/M^) to calculate it, while if it is given different from zero, Mw is not used by 
the program. Up to now no radiative corrections are included into the program. 

- The beam energy. 

- The degree of longitudinal polarization of the electron beam. 

- The number of requested events in the unweighted events sample. 

- The minimal and maximal allowed energy of the final particles. The maximal allowed 
energy of the photon can be modified by the program (and notified in the output file) to 
^max ^ - Ef"^" - E^'"" < Eb{l- ml/El) if the given maximum is too high to be 
consistent with the other cuts supplied by the user. 

- The minimal and maximal polar angles (in radians) of the final electron (9- = 9i) and 
positron = tt — $2) with their initial directions. 

- The minimal angle between final lepton and photon (it can be set to zero). 

- The minimal angle between initial lepton and final photon (it can be set to zero). 

- The maximal allowed acoUinearity (acoUinearity = \ tt — 9i — $2 |, where 9i and ^2 are 
the angles between final electron and positron momenta with the initial direction of the 
electron) and acoplanarity (acoplanarity = | tt — (02 — (f>i) |), 0i and 02 are the azimuthal 
angles of the final electron and positron. 

- A flag which allows writing (1) or no writing (0) the generated events into the created 
file BHAGEN-IPH.EVN. 

- A flag which allows writing (1) or no writing (0) the separate channel contributions. 

- A flag which allows to take into account separate contributions from s- and t-channel, 
which can be useful for tests, but require the switch to be set to 0, for comparing the results 
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with experiment. If the flag is equal to fuU version of the program is used; if it is equal 

to 1 only t-channel photon exchange contribution is included (with up-down interference 
neglected); if it is equal to 2 only s-channel contribution is included (from both 7 and 
exchange diagrams). 

- A seed to initialize the pseudo random number generator. 

Two output files are created. The first one BHAGEN-IPH.RES contains the calcu- 
lated cross-sections for both weighted and unweighted event sample and if requested the 
contributions from all ten channels. The second one BHAGEN-IPH.EVN contains the 
generated events (if the proper fiag was set to 1). Both files report the parameters, which 
are read from the input file. 
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TEST RUN INPUT. 



Input file to run with VMS system: 



$ruii bhagen-lph.exe 
91.1887 
2.49661 
80.304 
. 224482 
47.5 
1.0 
10000 
4.75 
23.75 
23.75 
0.7 
0.8 
0.01 
0.1 
180. 
1 





Z mass (GeV) 
Z width (GeV) 

W mass (GeV) (used only if sin**2(thW)=0. ) 
sin**2(t]iW) (if=0 then it is calculated =1-(M_W/M_Z)**2 
e_beam (GeV) 

longitudinal polcirization of the initial electron 
nevent - numb, of events in unweighted event sample 



23.75 
47.5 
47.5 

2.3 

2.2 



egjnin , eg_max (GeV) 
el_min, el_max (GeV) 
e2_min , e2_max (GeV) 
thl_min,thl_max (rad) 
th2_min,th2jiiax (rad) 
thlg_min=th2g_min (rad) 
thg_min (rad) 

acolcut (deg) ,acoplcut (deg) see comment (8) 
printing(l) no printing(O) generated events 
printing(l)/no-printing(0) info on separete channel contrib. 
s-t channel switch see comment (9) 

seed to the random number generator 



see comment 
see comment 
see comment 
see comment 
see comment 
see comment 
see comment 



(1) 
(2) 
(3) 
(4) 
(5) 
(6) 
(7) 



180. 



123456789 
$exit 
Comments : 

(1) - min. and max. photon energy in GeV. The minimum photon energy 
cannot be set to ! ! ! ! The maximum photon energy is set to 

eg_max = e_beam * (1 - me**2/e_beam**2) in case it is given bigger 
then that value. This is the biggest value of the photon energy 
allowed by four momentum conservation. 

(2) - min. and max. final electron energy in GeV. In case el_min=0. 
it is set to the electron mass by the program. 

(3) - min. and max. final positron energy in GeV. In case e2_min=0. 
it is set to the electron mass by the program. 

(4) - min. and max. electron polar angle measured from initial 
electron direction. 

A user is kindly requested not to use this program with 
thl_min < 0.001 rad. as the formulae used are not adequate 
for that region. 

(5) - min. and max. positron polar angle measured from initial 
positron direction. 

A user is kindly requested not to use this program with 
th2_min < 0.001 rad. as the formulae used are not adequate 
for that region. 
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(6) - minimal angle between final electron (positron) and the photon. 
In case they are not equal the smaller should be given. 

It can be set to 0. 

(7) - minimal angle between initial electron(positron) and the photon. 
In case they are not equal the smaller should be given. 

It can be set to 0. 

(8) - maximal acolinearity and acoplanarity allowed. 

Where acolinearity is defined as I (pi - thl - th2) I , and thl(th2) 
is the angle between final electron (positron) and initial electron. 
While acoplanarity is defined as I (phi2-phil-pi) I , where 
phil and phi2 are respectively the electron and positron 
azimuthal angles . 

(9) - 0:all included; 

l:t-chcinnel QED only with up-down interference neglected 
suitable for tests at small angles 

2:s-channel only, both photon and Z boson exchange include 



Additional, more detailed, event selection can be done by modification 
of the subroutine TRIGGER in the FORTRAN code. 

The cross section given in the output corresponds to the complete event 

selection e.g. with the cuts imposed by both the input parameters 

in this file and the cuts imposed in the subroutine TRIGGER. 

A simple example how to use this subroutine is given below. 

An event with electron having energy between . 1 GeV and 10 GeV is 

rejected and accepted if the electron energy is outside this region. 

c 

c This is a trigger which is to be defined by user allowing for the 

c rejection (k=0 has to be set) or acceptance (k=l has to be set) 

c of the generated event . 

c The final particles energies and angles are stored in the 'ev' vector. 

c They are given in CM frame of initial electron and positron, with 

c the z-axis along initial electron momentum. 

c 

c 

c ev(l) - electron energy in GeV 

c ev(2) - theta (electron) in rad. 

c ev(3) - phi (electron) in rad. 

c ev(4) - positron energy in GeV 

c ev(5) - theta (positron) in rad. 

c ev(6) - phi (positron) in rad. 

c ev(7) - photon energy in GeV 

c ev(8) - theta (photon) in rad. 

c ev(9) - phi (photon) in rad. 
c 
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subroutine trigger (ev,k) 

c 

implicit real*8 (a-h,o-z) 

dimensioii ev(9) 

c 

if ((ev(l) .gt.O.ldO) .and. (ev(l) . It . 10 .dO) )then 

k=0 

else 

k=l 

endif 

c 

return 
end 
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TEST RUN OUTPUT. 



Output file BHAGNE-IPH.RES. 



1.0000 
47.5000D+00 GeV 
47.5000D+00 GeV 
2.37500D+01 GeV 

< 2.30000D+00 

< 2.34159D+00 



BHAGEN-IPH, VERSION OF 22-JAN-1996 
Z mass (GeV) : 91 . 188700D+00 
W mass (GeV) : 80 . 304000D+00 
sin**2 (th_w) = 0.224482000000000 
Total Z width (GeV) = 2.4966100D+00 
E_beam = 47.5000 GeV 
Longitudinal polarization of electron 
23.7500D+00 GeV <= electron energy <= 
23.7500D+00 GeV <= positron energy <= 
4.75000D+00 GeV < photon energy < 

7.00000D-01 rad < theta (electron) < 2.30000D+00 rad 
9.41593D-01 rad < theta (positron) < 2.34159D+00 rad 
Minimal allowed angle between photon and final lepton = l.OOOOOD-02 rad. 
Minimal allowed angle between photon and initial lepton= l.OOOOOD-01 rad. 
acollinearity < 3.14159D+00 rad. 
acoplanarity < 3.14159D+00 rad. 

no. of events with weights > max. weight: 

Total cross section of the events with weights > max. weight : 

cross section (nb) = . OOOOOOOD+00 +- O.OOD+00 

Number of generated (unweighted) events in the sample: 10000 
Number of hits (weighted events) in the sample: 130202 
Summary : 

Total cross section obtained using weighted events : 
total cross sectionl (nb) = 1 . 3745087D-02 +- 3.61D-05 
Total cross section obtained using unweighted events: 
total cross section2 (nb) = 1 . 3874919D-02 +- 1.39D-04 
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Output file BHAGNE-IPH.EVN - only first two events are reported. 



BHAGEN-IPH, VERSION OF 22-JAN-1996 
Z mass (GeV) : 91 . 188700D+00 
W mass (GeV) : 803.04000D-01 
sin**2 (th_w) = 0.224482000000000 
Total Z width (GeV) = 2.4966100D+00 
E_beain = 47.5000 GeV 
Longitudinal polarization of electron 
23.7500D+00 GeV <= electron energy <= 
23.7500D+00 GeV <= positron energy <= 
4.75000D+00 GeV < photon energy < 

7.00000D-01 rad < theta (electron) < 2.30000D+00 rad 
9.41593D-01 rad < theta (positron) < 2.34159D+00 rad 
Minimal allowed angle between photon and final lepton = l.OOOOOD-02 rad. 
Minimal allowed angle between photon and initial lepton= l.OOOOOD-01 rad. 
acollinearity < 3.14159D+00 rad. 
acoplanarity < 3.14159D+00 rad. 
Structure of the event : 



1 . 0000 
47.5000D+00 GeV 

47.5000D+00 GeV 
2.37500D+01 GeV 

< 2.30000D+00 

< 2.34159D+00 



Electron energy (GeV) 
Positron energy (GeV) 
Photon energy (GeV) 



3.83488D+01 
4.26590D+01 
1.39922D+01 
4.35498D+01 
4.54257D+01 
6 . 02442D+00 



1 . 33469D+00 
2.13715D+00 
1.02513D-01 
1.14156D+00 
2.12109D+00 
3.63628D-01 



theta(electron) (rad.) 
theta(positron) (rad.) 
theta(photon) (rad.) 
1.53144D+00 
4.65588D+00 
5.11897D+00 
8.91211D-01 
3 . 98290D+00 
5 . 15577D+00 



phi (electron) 
phi (positron) 
; phi (photon) 



(rad.) 
(rad.) 
(rad.) 
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